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AISTRACT 


Vm hav* atatlatlcally analyiad tha dlaparaton naasuraa of a aampla of U9 
pulaara In tha tnnar Galaxy (|t| < 50^) to daduca tha larga-acala dlatrlbutton 
of fraa tharaal alactrona tn thla ragton< Tha dlaparaton aaaaure dlatrlbutton 
of thaaa pulaara ahowa atgnlflcant avldanca for a dacraaae In tha alactron 
acala halght froa a local valua greatar than tha pulaar acala height to a 
value lean than tha pulaar scale hef.^ht at galactocantrtc radii Inalde of 7 
kpc. An Increase In tha alactron dantlty (to a valua around •IS cn~^ at 4*>5 
kpc) nuet accoiapany such a dacraaae In acala hatght. Thera la also evidence 
for a larga-acala warp In the alactron distribution below the b ■ 0° plane 
Inalde tha Solar circle. Wo propose a nodal for tha electron distribution 
which Inci’poratas these festuras and present Monte Carlo ganarsted dispersion 
measure distributions for parametara which bast reproduce tha observed pulsar 


distributions. 
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PRECEDING PAGE BUNK NOT FILMED 

1. INTRODUCTION 

Pulsar radio signals ara dlsparsed by froa alaetrons In tha Intarstellar 
medium which lie along the line of sight. Pulsar dispersion measures are 
therefore good probes of the Interstellar electron density and Its galactic 
distribution. Unfortunately, few pulsars have reliable Independent distance 
determinations and. since the dispersion depends both on distance and c^lectron 
density, the mean electron densities along lines of sight to Individual 
pulsars cannot be determined In most cases. However, statistical 
distributions of pulsar dispersion measures can be used to study the electron 
distribution, provided something Is known or assumed about the distribution of 
pulsars In the Galaxy. 

Traditionally, the opposite approach has been taken. The pulsar galactic 
distribution Is determined by assuming a simple model for the free electron 
distribution which la in agreement with other observations (Davies, Lyne and 
Selrsdakls 1977, Taylor and Manchester 1977, Lyne 1981). These other 
observations are free-free absorption of the non-thermal radio continuum. 
Interstellar scattering of extragalactlc sources, HI absorption of pulsar 
signals, and radio recombination lines. Both free-fvee absorption of 
galactic radio vr»ectra at low frequencies and Interstellar scattering Indicate 
a Z scale height for the absorbing electrons of 500-1000 pc. (Bridle and 
Venugopal 1969, Readhead and Duffett-Smlth 1975). The best estimates of mean 
electron densities on klloparsec scales come from the distance determinations 
to pulsars usi>.g HI absorption (Abies and Manchester 1976; Helsberg, Rankin 
and Borlakoff 1980). From the distances and dispersion measures for a sample 
of low latitude pulsars, Welsberg st al. (1980) conclude that the mean 





olpctron (i^nalty In the pl«n« 1* <ng> ■“ .025-. 03 c«”^ In dlrectlone away from 

the galactic center. Therefore, the "etandard" model haa been an exponentl. 1 

or Gauaalan Z dependence with ecale height ~ 1000 pc and central denalty, n^ 

-1 

.03 oil . However, the HI abaorptlon meaaurementa alao give mean daneteleH 
that are aomewhat higher toward the Inner Galaxy, within a longitude range 
130 1 ^ 30”, auggeatlng that n^ te dependent on galactic redlue. Radio 

recombination line obeervatlone (Lorkman 19/6) alao Indicate, from the radial 
dletrll.it Ion of lonleed gaa, that electron denaltlea Increaae In t^'a Inner 
Calaxv. There la evidence then, that more complex modela may he needed to 
dcacrthe rite galactic dlatrlbutlon of electrona. 

In ihia papi r, we aiualyce *'he pulaar dlaperslon meaaure diatrlbutlone In 
the lnnt‘r Galaxy to atiidy the variation of the denalty and acnle heighc of the 
elecironH with g.ilac rocent rlc radlua. Wc uae reeulta of provloue analynea to 
make Haaiimptlona about the pulaar galactic dlatrlbutlon. Modela for the 
electron dlatrlbutlon can thee be tested by generating Monte Carlo dlaperslon 
measure distributions and comparing these to the obaerved distributions. We 
show that a particular density and scale height variation with R and Z will 
repr >duce the observed pulaar dispersion measure dlatrlbutlons, and dtscusa 
Mic consistency of this type of model with other observations of free 
Interstellar «»lev'tr ms. 


II. PULSAR niSPERSION MEASURES 

The dispersion measure (DM) la defined as the Integral of the electron 
density along the line of sight to a pulsar; 

“ r* n (s) ds 

o c 


DM 


( 1 ) 
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It is dctaralnad fro« the tnturvala batwaan pulaa arrival cluaa at diffarant 
fraquanctaa and la tharafora a dlract obaarvattonal paraaatar. Inforaaeton on 
eha Z dlatrlbuCion of alaetrona and pulaara can ba obtalnad by plotting the 
nuabar of pulaara In Intarvala of A ■ OH alnb, tha Z coaponant of dlaparalon 
aaaaura. Tha obaarvad dlatrlbutlnn can ba coaparad to tha axpactad 
dlatrlbutlon of N(A) varaua A darlvad by aaaualng that tha pulaara «nd 
alaetrona are funettona of Z only* If for axaapla both have axponantlal 
dlatrlbutlona In Z, 


n (Z) - n axp (-Izl/h ) (2) 

e o e 

N (Z) - N exp (-|Z|/h ), 

P o P 


than the axpactad form of N(A) la: 


N(A) 


N , TT 

n nr ^ 

o o e 

0 


-1 


• A < "o^ 


, A > n h 
o e 


(3) 


(Taylor and Nancheater 1977). This N(A) dlatrlbutlon will have three dlatlnct 
shapes, depending on t' relative sizes of the electron and pulsar scale 
heights, hg and hp. If (1) h^ > hp, then N(|A|) la aonotonlcally decreasing 
out to |A| - n^jhg (2) If hg “ hp, then N (|A|) la constant out to |A| - n^jh^ 
(3) if hg < hp, then N(|A|) la monotonlcally Increasing out to iA| " n^h^ 
where It abruptly drops off. Tha general dependence of the shape of N(A) on 
the ratio h*/h_ was noted previously by Gould (1971) and by Terzlan and 

B p 

Davidson (1976) and does not change with different functional forms for the 



t<!leccron and pulaar dUtrtbutlona, •« long aa thay ara both aonotonteally 
decreiialng fiinctlona of 7. • 
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A hlatograa of N(|A|) varaua |A| for all tha known pulaara tuma out to 

be a decreasing function of |A|, with “ ” 0^0 ~ P** » Indicating 

that h^^, hp throughout moat of tha Calaxy (Taylor and Manchaatar 1977, 

llardliiM '®8l). Equation (9) can ha fit to tha obaarvad dlatrlbutlon to 

determine one of the quantities h^, h^ or n^, aaoumlng values for tha other 

two. This approach has been used for the Molonglo survey (Manchester 1979) to 

-3 

give hp “ 390 pc. for assumed values of h^ ■“ 1000 pc. and n,j “ .03 cm 

Oallly et al. (1978), u (ng distributions for electrons and pulsars which were 

Intermediate between Cat islan and exponential, obtained ^IZp|> “ 380 pc. for 

7 SO < vl7.„h < 1900 pc. and n^j - .029 cin"^. 

Xlu***** It? cvldtMice, however, that the electron Z distribution Is not so 

simple. Komosarolf et al. (1973) noted an asymnatry In a plot of IW versus b. 

In which pulsars .at southern latitudes had higher dispersion measures. They 

suggested that a model of n^{Z) consisting of a high density electron disk 

with n^, - .12 em"^ with s small thickness of ' 124 pc. would explain the 

effect if the Sun la displaced 22 pc north of the center of the disk. Harding 

(1981) (ound a similar anymmetry In a plot of N(A) versus A and proposed a two 

-3 

component model consisting of a thin disk of electrons with n^ " .10 cm of 
half thickness AZ - 40 pc. and a large scale height exponential component with 
n^, - .03 cra"^ and h^, - 1000 pc. If the Sun were displaced 20 pc. to the north 
of this distribution, then equal pulsar scale heights of 340 pc. were obtained 
for both -t-Z and -Z pulsars . rora fits to the N(A) versus A histogram. It was 
also noted ( his asymmetry In the N(A) versus A distribution was moat 

pronounced toward the galactic center, e.g., for a sample of pulsars 
restricted to a longitude range |l| < 90° (310° < I < J0°). This suggested a 
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poasibl* radial dapandanea of tha alaetron dlak danalty If eha dlak la 
raaponalbla for ■oat of tha aayaaatryt 

Wa now eonatdar tha diaparalon aaaauro dlatrlbution of thla aaapla of 
pulaara with 111 < SO** in aora dotal !• Plgura 1 ahowa a aeattar plot of ON 
olnb ( ■ A) (vortical DM) voroua DM eocb (horicontal DM) for thla aaiaplo of 

149 puloaro. Thla plot ravaala that aaaan^.lally all of tha aayaaetry In tha 

s .<1 

N(A) varaua A dlatrlbution la dua to pulaara with DM eoab - 100 pe« ca 
l.a., aora dlatant pulaara. Tha DM slnb dlatrlbution of naarby pulaara la 
actually qulto ayaaatrlc. Tha aaan valua of DM ainb appaara to dacroaaa with 
DM coab and roach a alnlaua at DM coab ^ 200 pc. ca~^. Tha dlatrlbution aloo 
narrowa with Incraaalng DM coab. 

In Pigura 2, the aaapla la dlvldad Into groupa having DM coab > 100 
(dlatant pulaara) and DM coab < 100 (nearby pulaara) and plottad In the 
hlatograa fora N(A) veraua A, for the entire Inner Galaxy, lit < 50°, and for 
the central Inner Galaxy, |i| < 30**. To sake a unlforo aaaple, only pulaara 
detected by the aecond Molonglo aurvey (Mancheater et el. 1978) have been 
plotted In Plgure 2. The more aenaltlve Areclbo aurvay (Hulae and Taylor 
1974) covered a aoall longitude range, 42** < t < 60**, and the pulaara detected 
only by thla aurvay have been removed from the aaaple. 

In aplte of the atatlatlea, a algnlflcant difference In ahape, width and 
mean la apparent for dlatrlbutlona of nearby ani dlatant pulaara. The |i| < 
50° hlatograa for nearby pulaara la alallar In ahape, width and aaan to the 
dlatrlbution of pulaara In the outer Galaxy (|t| > 50**). Within atatlatlcal 
errora. It la ayaaetrlc and centrally peaked with A^^j^ = n^jh^ ~ 25 pc. ca"^. 
According to the model of equation (2), <h^> > <t.p> for thla aaaple of 
pulaara. If <h^> ■ 800 pc., then <n^j> - A^x^^e^ ^ which la 

conaiatent with aodela which have been derived froa the entire aaaple of known 
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piilsorn. The Ut < 50^ hiatogriw for dlecant puleara la docldadly aajTMattic 
with « tiean - - 4 pc. cm~^, and a ahapa laora auggaative of tha eaaa <h^> 
<hp>, where N(A) la conatant or Incraaalng out to dgm,|« The hlatograiaa for 
Id < 30** are laaa centrally peaked than the 111 < SO^ hlatograaa. Bapeelally 
arrlklng la the doubly peaked dlatrlbutlon for the aaaple having DH coab > 100 
with a mean also around -4 pc. cm . The dlatrlbutlon la fairly aymmetrlc 
/‘bout this moan and tha shape auggeata that <hg> < <hp> for thla aanple of 
pulaafH, with - A^l ~ 15 pc. If <hg> 200 pc., algnlf icantly 

less than a pulsar scale height of 340 pc., than <ng> '• .08 cm . Tharefore, 
to account for the shapes and widths of the obaerved dlatrlbutlons, the 
electron scale height must decrease significantly and the density must 
Increase toward the galactic center. 

The displacement of -4 pc. cm In the distant pulsar distributions may 
be caused by a dtsplacemont In the electron or pulsar distributions with 
respect to the plane of the Sun or, most likely, both. Since the asymmetry 
does not appear strongly for the nearby pulsars, the displacement must be much 
greater In the inner Galaxy than it Is locally. This explanation Is 
consistent with a warping of the galactic plane toward negative Z, an effect 
which has been seen In CO (Cohen and Thaddeus 1977, Solomon, Sanders and 
Seovllle 1979), 'III (Loc.kman 1977, 1979) and HI (Quiroga 1974) data. From 
tliese observations, the maximum displacement south of the plane seems to occur 
between galactic radii of 6 and 7 kpc and 1 e about 40 pc. in magnitude. This 
seems roughly consistent with the displacement needed to account for a DMslnb 
difference of 4 pc cm , since If <ng> '« .08 cra”^ then AZ ~ Ag/<ng> 50 pc. 

From the simple arguments outlined above, a free electron distribution 
that Is consistent with the obaerved distributions of oulsar dispersion 
measures In the Inner Galaxy should have the following properties: 
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1) An •ffaetlvn «c«l« hnltht iihlch 1« • function of t, Inerooclng from « 
value laao than tho pulaar acala halght to a local value graator than the 
pulaar acala height 

2) A danalty trhlch la a dacraaslng function of R, at laaat batwaan 6 and 
10 kpc. 

3) A Z dlatrlbutlon whoaa naan la nagatlva for R i 8 kpc. 

III. KODBL FOR THE RLBCTRON DISTf TBUTION 

We have teatad k'paelfle nodaia for the electron dlatrlbutlon In the Inner 
Galaxy by neana of a Monte Carlo progr*^ which calculatea dlapcralon meaaurea 
to randomly choaen polnta and ganerataa N(A) veraua A hlatograaa to coapare 
with the data. The polnta are choaen fron an aaauaad pulaar dlatrlbutlon in 
the Galaxy and obaervatlonal aelactlon ef facta are taken Into account. A 
particular functional form wae chop**n for the electron Z dlatrlbutlon and tne 
paraaetera were varied to obtain disperalon aeaaure dlatrlbutlona which were 
C'lnalatent with the obaerved dlatrlbutlona. 

a) The Pulaar Distribution 

The program populates the Galaxy w'^th pulaara ualng a Monte Carlo 
method. In galactocentrlc, cylindrical coordlnatea the radlua R, Z distance, 
and angle, 6, for each pulsar are chosen according to specific probability 
functions. These functions are chosen to be consistent with determinations of 
the pulsar galactic dlatrlbutlon (Lyne 1981, HarJlng 1981). The probability 
of choosing a pulsar as a function of 6 Is taken to be constant, l.e., 

<iP0 

■ constant. Because we assume aslmuthal symmetry about the galactic 
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center for both the puleer and electron dlatrlbutlona and baeauae alaoat all 
pulaara on the oppoalta aida of the Galaxy froa '’ha aolar aystaa ara 
unobeervable, we need to populate only one quadrant of tha Galaxy with 
pulsars. Hence, 6 Is allowed to vary only between 0 and «/2. 

The probability of chooslnf a pulsar aa a function of R le taken to be 

- R exp (-IR^-Rl/W) . (4) 


W and Rq are chosen Co be 3 kpc and 5.5 kpc, respectively, such chat the peak 
pulsar surface density projected onto Che plane occurr at a radius equal Co 
5.5 kpc and Is 4 times Che surface density at Che location of the solar 
system, which is assumed to be at R ■ 10 kpc. This probability function is 
integrated from R > 0 to R ■ 10 kpc and a correspondence table la set up 
correlating 's’ues of the Integrated probability with valuea of radius. 

The probability of choosing a pulsar aa a function of Z is taken to be 



^ exp (-|Z|/hp) 


(5) 


where Is Che pulsar scale height. A constant pulsar scalci height of 340 
pc. has been chosen for nearly every toudel Cested. 

The probability of choosing a pulsar aa a function of radio luminosity is 



( 6 ) 


where Is the lower luminosity limit. This choice of probability 
distribution corresponds Co recent, calculated pulsar lumlnoslC/ functions 
(Manchester 1979, Lyne 1981). The derived lower luminosity limit is around .3 
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■Jy kpc^i howavnr, • larger value of • 1 mJy kpe^ vaa uaed whan teatlng 
(Bodala. Thla graatljr apaada up progran asoeutioni and pulaara wish 
luainoatetaa lane Chan 1 «Jy kpe^ have only boon doCaetad at dlaCaneaa Iooh 
chan ' 300 pc> In the Molonglo aurvay* The raaaon for ehooaing a lualnoaicy 
for each pulaar la Co chock within cha prograa for doCacCablllty by the 
Molonglo aurvay. Pulaara noc obaarvabla by thla aurvay are dlacardad by the 
prograa. 

b) The Electron Dlacrlbutlon 

The nodal we conalder for the electron dlotrlbutlon conalate prlnarlly of 
CtM coaponenta: the dlak alectrona and the extended elictrona. It la 

deacrlbed by the following functional fom: 


n®(R), |Z-Z^(R)I<AZ 

n (R.Z) - (7) 

* n^(R) expl-|Z-Z^(R)|/h^(R>), |Z-Z^(h>l>AZ 

Th? dlak coaponenC conalaca of a relatively thin layer of denaicy n^(R) lying 
In Che plane defined by Zg(R). The chlekneaa of the dlak coaponenc, AZ, was 
Caken to be Independent of radlua. Both the thlckneaa of the dlak and the 
dlak electron denalty, which waa taken to be a function of radlua, were varied 
fron one trial to the next to aearch for a viable nodal. The electron denelty 
within the dlak waa held conatant with Z. 

The extended component conalata of free electrona extending above and 
below the dlak oonponent. The denalty of the extended conponent was taken to 
fall off exponentially with dlatance froa the galactic plans. Thua, the 
extended conponent la characterized by a central denalty, n^CR), and a acale 
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height, h,(R). at each galactocantrtc radlua. Dlffarant functional foraa for 
both tha acala halght and tha eantral danalty wara taatad by tha prograa> 

Both coaponanta hava an R dapandant aaan, Z^(R)t irtilch dafinaa a warp In 
tha plane of tha alactron dlatrlbutlon, and which la choa<m to hava tha 
functional fora: 

Z (R) - -Z expC-(R-R )^/AR^J (8) 

o' max ' o ' 

We hava adopted the valuee - 6.5 kpc. and AR " 1.7 kpc., to be conalatent 
with the warp aean in HII reglona in the Inner Galaxy (Lockman 1979). 
c) Tearing Procedure 

I1te cylindrical coordlnatea and lumlnoalty of a saaple pulaar are chosen 
according to the probability distributions described In Ilia. The longitude, 
latitude and distance d of the sample pulsar to the Sun are then calculated. 
With an assumed electron distribution, the dispersion measure la determined by 
Integrating along the line of sight from the pulsar to the Sun. The flux 
density, S, is determined from the luminosity and distance by S •• Ld end Is 
compared to the minimum detectable flux density at that 1, b, and DM for the 
Molonglo survey. The effects of a variable sky temperature from the radio 
continuum background and the dispersion measure cutoff of the survey were 
taken Into account in determining the minimum detectable flux density (see 
Taylor and Manchester 1977). The latitude, longitude and dispersion measure 
of each detectable pulsar are accumulated In a file In order to build up a 
statist leal sample which reflects the model being tested. 

We can then plot various distributions of N(A) from the Monte Carlo 
generated sample and compare them in shape, width and mean with the observed 
distributions. Different variations of n^(R), n^(R), h^(R) and 
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■quaclon (7) iftrc taaead until a nodal waa found which t^va a'StlafacCory 
atraanant with tha data. 


IV. RESULTS 


A nodal was found which raproducaa the data In Figure 2 raaaonably 
wall. Tha N(A) dlstrll<utlona producad by thla nodal era shown In Figure 3. 
Tha functional foms for tha variation of alactron danaity and elactron scale 
height with galactic radius ara shown iu Table 1 and plotted in Figure 4. 

It Is evident that the naln feature of the nodel Is the drsnatlcally 
decreasing electron scale height toward the galactic canter. In order that 
the model produce dispersion measure distributions having tha characteristics 
of the observed distributions with DM cosb > 100, it is essential that the 
electron scale height be less than the pulsar scale height at galactic radii 
less than about 7 kpc. The radius at which the transition h^ > hp to h^ < hp 
occurs Is ■ 7.5 kpc and Is a critical paramb^sc of the nodel. It may also be 
Important for the central density of the extended component to Increase 
rapidly with decreasing R In the vicinity of this transition region. In the 
model described In Table 1, this central density then levels off and remains 
relatively constant with decreasing R Inside about R " 7 kpc. It may be that 
this sudden Increase In electron density corresponds to the location of a 
spiral ariit, or to the beginning of a "galactic ridge" of electrons (Seacord 
and Gottesman 1977, Lockman 1976, Gordan and Gottesman 1971) seen in radio 
recombination line surveys. Models without this rapid Increase In electron 
density (e.g., a linear Increase from R > 10 kpc to 4 kpc) required a linear 
decrease in electron scale height from 1000 pc. to x 200 pc. over radii of 10 
kpc to 7.5 kpc. These models produced >rer fits to the N(A) histograms for 


luiMrby pulii«rii bacaus* Cha avaraga aeala halght batwian 10 and 8 kpo waa too 
amall* 

Tha warping of tha alactron dlatrlbutlona producaa tha horlaontal offaata 
In Che N(A) hlstugrama for DM eoab > 100. In Cha nodal which produead cha 
hlacograins In Figure 3, Cha maxlmun dlaplacananC • balow Cha p1«na aC 6.3 

kpc (cf. Squaclon (8)) waa 120 >c. A of 60 pc. can produce Cha aame 
offoec of '> -4 pc cit~^ If Cha Sun la dlaplacad 20 pc. above Cha 7. • 0 plana of 
Che electron dlatrlbutlon. Thla dlaCorclon la more in agraananC with Cha 
Inner Galaxy warp eeen In ocher gaaaoua Craeara. Furchamora, a naxlmun dla- 
plao.eracMtc aa large aa 120 pc. aC 6.5 kpc. nay be InconpaClbla with recon- 
blnaclon lino data in the Inner Galaxy (Locknan, private communication). 

Thore In no iinauranco tl»at the model In Table I (a unique. Ocher model n 
with Homewhac different comblnitloua of h.!', «.sd N^(R) night alao reproduce 
Che data. The fit to cho data waa very aanalclva to certain of the paranetera 
while for other paranetera the valuea choaan made relatively little 
difference. It la felt that Che choice of paranetera for galactocenCrlc radii 
loan tliun 4 kpc. la not particularly inportant aa relatively few ;nilaara In 
the date or delectable polnCa In the model sample thla region of the Galaxy. 

An .1 rcniilt, eloctron denaltlen and Hi'ale heights are poorly determined by tlie 
model for thin region. The model la also relatively Inaenaltive to tlie 
electron disk denalty Inside about 6 or 7 kpc. The lines of sight to only a 
few sample pulsars pass through tills portion of Che disk, and hence only a 
small percentage of dispersion measures are affected by the inner disk. The 
two ooinponent form adopted for the electron Z-dlstrlhatlon |cf. Rquatton (7)| 
was one of nkit homaC leal convenience and computational efficiency. Any other 
set of funi'tlons giving a similar profile in Z, for example a Gaussian disk 
plus .in exponential extended component, would have worked an well. 


IS 


Mod* If eonfUclng of nn fxCfndfd eoaponont Which !■ tndai^ndent of R tforc 
daflnlCtly unfuccfsaful In roproducl'ts data, avan vlth R dapandant dlak 
eoaponanta. By way of llluatratlon, wa ahob tha N(A) dlatrlbutlona raaultlng 
froa cha following alactron ac>dals 


" (^*R) " *025 exp(-*|z|/h ) -f nf (R) (9) 

V O V 

where 

.0l(l+l.5(10-R)l , Z<AZ 

n“(R) - 

• 0 , Z>AZ 

and AZ " 60 pc., h “ 1000 pc. Tha dlak danalty In thla nodal haa a linear R 
variation which la very alallar to that In Figure 4, but tha extended 
exponential component haa a central denalty and scale height which la constant 
with R. Thla model Is also similar to one recently adopted by Lyne (1981) to 
determine tha pulsar galactic distribution. Figure 5 shows the resulting N(A) 
distributions, which clearly have different shapes and wllths froa the 
observed distributions. Ue have not Included a warping of the electrons in 
this model, no that tha negative offset In the data has, of course, not been 
reproduced. 

As an approximate test of the relative "gocdness or fit” of the 

distributions generated by different electron models with the data, we have 
2 

determined tha x values and probabilities for each computed histogram. These 

are listed In Table 2 for the best fit modal, the constant scale height model, 

2 

and the simple exponential model. By tha x teat criterion, the best fit 
model provides a significantly better fit than tha simpler models, especially 
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for chr |t| < 30**, DHcoiib > 100 aaaplo of pulaara. Boetua* te la difficult to 
evaluate the lumber of Independent paraaetera In a Model aa eoaplax aa oura, 
the nuMber of depreea of freedoa are aoaawhat uncertain ao thla type of teat 
cannot be conaldered definitive. We therefore have not atteapted to uae 
HtHtlatlc.il tentH to put uncertalntioH or any of the aodel paraaetera. 

However, we have ahown here that the obaerved dlatrlbutlona, daaplta the email 
minberH of pulaara, differ aignlf Icantly from thoae derived from the almpler 
modele (l.o., the probability that the true parent dlatrlbutlona could be 
those p,eneratod by the almpler modela la quite amall). 

Another obrervatlonal conatralnt by which electron aodela can he tested 
for cons'lHtency with the data Is the distribution of number of pulsars versus 
DM cosh. The ohmirved distribution Is a decreasing function of DM coab, 
reflecrlng the sensitivity limit on detection of distant pulsars. Provided 
that one c.^n ad«-quately account for selection effects, a viable electron model 
should give a ratio of detectable pulsars with DM cosb < 100 to those with DM 
cosb > 100 that Is In agreement with the data. Since the Molonglo Survey was 
sensitivity limited rather than dispersion measure Halted (the maximum 
dtitectahlo DM was 780 pc cm In the plane) and there are only a few pulsars 
In the nb«erved sample with DM > 400 pc cm , a viable model should also not 
produce too many very high dispersion pulsars. Our best fit model glvea 
ratios of N(DM cosh < 100)/N(DM cosb > 100) which are In good agreement with 
the data for both Ul < 30® and |t| < 50®. The model also produces a sample 
with '"5.5 percent having DM > 400, In comparison to the observed sample with 
^ 3 percent having DM > 400. By contrast, the model of Equation (9) giving 
the dlstrlbut lens shown in Figure 5 produced about twlca the ratio N(DM coab 
< lOO/iNfOM cosh > 100) «s la seen In the data. Since the pulsar luminosity 
cutoff, Lq, Is probably less than the value of 1 mJy kpc we have assumed, and 
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« lower Lq would glvo an even larger ratio of naarliy to dtatant pulaara, thl* 
■odal aaaaa to be quite Ineonalatant with the data. 

Since the focua of thla work waa to uaa the pulaar dlaparalon neaaure 
data to Inveatlgata the electron danalty dlatrlbutlon, wa have bald the pulaar 
paraaatera In the nodal constant. Our beat fit model la relatively 
Inaen'iltlve to changaa In the ratio of pulsar aurfaca densities at 5 and 10 
kpc, at least with the exponential form wa have assumed. The low luminosity 
cutoff In the pulsar luminosity function. L^, does not affect the distant 
pulsar distributions at all, since pulsars with low luminosities are only 
detectable within a kpc or so (depending, of course, on the value of L^). It 
does affect the distributions of nearby pulsars. In that low luminosity 
pulsars tend to fill In the central part of the N(A) distribution since they 
are not fully sampled In Z. For this reason, we have concentrated primarily 
on fitting the distributions for DM cosb > 100, and only the general outer 
shapes and widths of the distributions for DM cosb < 100. 

One could probably contrive a model where drastic changes In the pulsar Z 
distribution In the Inner Galaxy would entirely account for the observed 
changes In the dispersion measure distributions, without Invoking any radial 
dependence In the electron distribution. This model would require that 
pulsars be somehow excluded from the central galactic plane and peak In 
density both above and below the plane Inside of 8 kpc. We consider this 
solution to be physically unrealistic, as no Population I objects have been 
observed to exhibit this structure. 

It Is possible, though, that the pulsar scale height Is a function of 
galactocentrlc radius. However, the pulsar scale height cannot decrease very 
much toward the galactic center. As we have seen, the electron scale height 
must be less than the pulsar scale height over part of the Inner Galaxy. As 
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Che eleccroit BCiil* height decreaeea, Che electron denelcy auac tncreeae to 
produce the obaerved eaoiinC of OH alnb. But ON coeb would becoae Coo large 
for a significant fraction of pulsars ((hose In the inner Galaxy lying clone 
to the plane) If the electron density irare Increased too «uch> 


V. DISCUSSION 


We have ahown evidence Chac Che dlsperulon SMasure distribution of 
distant pulsars In the Inner Galaxy differs significantly from Chat of nearby 
pulsars, Indicating that the electron Z distribution undergoes major 
qualitative changes between the Sun and the galactic center. We have 
presented a model for the radial dependence of Che electron distribution which 
la consistent with the pulsar data. Tt should be stressed Chat the particular 
model presented here la not meant Co be Che only consistent model nor Is It 
meant to be correct In every detail. Rather, It Is Intended Co Illustrate Che 
general characCerlsClcs of a class of models which are consistent with pulsar 
dlspersloi measures. 

The model we have proposed In this paper incorporates several new 
features not previously considered in electron density models. One of these 
Is an eloccron ocale height which increases with galactic radius inside Che 
Solar circle, from a value which Is smaller Chan the pulsar scale height. 
Previous models have assumed a constant electron scale height throughout the 
Galaxy. This model also includes a radial dependence of the electron lenslty, 
wl>lch is a necessary consequence of the scale height dependence. If the 
electron scale height decreases Coward the Inner Galaxy, then the density must 
Increase Co compensate for a loss of vertical dispersion measure which would 
not be consistent with the data. Finally, our model Includes a radially 



d«p«nd«nt dlspUcwMnt of tho oloctron piano towarda nagatlva Z, a faaturo 
which la nacaaaary to axplaln tha conoldorabla aaount of aayaaacry In eha 
vortical dloparalon MaBuro dlatrlbutlon of pulaara In tho Innor Galaxy* It 
la poaalblo that tho pulaar dlatrlbutlon haa a atallar faaturo, but It could 
not, by Itaolf, account for tho obaorved vortical dlaporalon neaaurc 
uayniiBetry. 

Sovaral of tha proportloa of thla oloctron dlatrlbutlon havo. In fact, 
boon auggoatod by othor typea of obaorvatlon or analyala. Cane (1977) found 
that a two conponent Z dlatrlbutlon of olectrona, a thin dlak plua an extonded 
component, waa needed to account for the free-froe abaorptlon of the low 
frequency radio background* The longitude dlatrlbutlon of mean electron 
denalty along the llnea of sight to pulaara with HI abaorptlon distances haa 
Indicated higher mean densities In the longitude range 30° < t < 330°. (Abies 
and Manchester 1976)* The radial dependence of the electron density In our 
model can account for this and gives mean denaltlea In these directions which 
are consistent with those observed* In analysing the dispersion measure 
distribution of the pulsars with Independent distances. Hall (1980) derived a 
mean density, <n^> » .048 cm“^, and an electron scale height, h^ ’ 264 pc. 

The disparity of these results with previous electron models can perhaps be 
understood ns Indicating that the electron scale height and denalty, as we 
have found, vary In different parts of the Galaxy. The Z asymmetry In the 
pulsar dispersion measure distribution has turned up In several previous 
analyses (Komesaroff et* al* 1973; Davies, Lyne and Selradakls 1977). The 
existence of the asymmetry In the Inner Galaxy r.nd its absence in the outer 
Galaxy (Harding 1981) was a preliminary Indication of a radially dependent 
displacement of the electron disk* 
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Having used asauaptlons about the galactic dlatrlbution of pulsars to 
derive properties of the electron distribution, «a should ask how the pulsar 
galactic distribution derived fro« this electron andal would differ frosi 
previous results. As dlscuased earlier, the results of the Monte Carlo 
calculation are not very sensitive to the fora of the pulsar distribution. Ue 
have not assumed anything about absolute pulsar densities, since only density 
ratloH entered :nto the calculation. Therefore, iinlesa the rederived ptilaar 
distributions change drastically with this electron model, the results wa have 
obtained will be valid. 

The major effect of this model would be to shrink the distance seals for 
pulsars In the Inner Galaxy. There would be little or no effect on pulsars 
within I kpc, since the local form of our electron distribution Is the sane as 
In previous models, but distant pulsars would be somewhat closer then 
prcvlounly thoup.lit. Since pulsar densities are almost entirely determined by 
the low luminosity pulsars lying within the nearest 500 pc. or so, only the 
form of the radial distribution %K>uld be expected to change, not the absolute 
densities. Distances to the very nearby pulsars, which determine densities 
and, ultimately, the pulsar birthrate In the Galaxy, depend critically on the 
local fluctuations In the electron density snd the presence of HIT regions. 

The lonl/.ed gas Is probably clumped on scales of I - 100 pc. (McKee and 
Osirlker l*)77, Dickey et al. 1981), hut this "fine structure' has been 
averaged over In our large-scale model. Local HIT regions affect somewhat the 
dispersion measure distributions of nearby pulsars, but make negligible 
contributions to the dispersion measures of more distant pulsars. 

In summary, pulsar dispersion measure) can provide a great deal of 
Infornatlon on the large-scale galactic distribution of Ionized gas. It la 
uncomplicated by a lack of knowledge of the temperature, fractional ionization 
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•nd B«all*aeal« atruecur* of eha gaa, aa ara othar probaa of tha alactron 

danalty. Tha only proaant obaarvatlonal datanilnatlono of <n^> com froa 

pulaar dlaparalon Maauroa; ochar obaarvaelona aaaaura <n^> froa ealaalon 

aeaaures. Coablnlng chaae two quancltiea can dataralna a vary iaportant 

paraaatar doaerlblng the atructura of the lonlaed gaai the filling 
2 2 

factor, f ■ <n > /<n >. Higher electron denaltlea In the Inner Calaxy, aa 

V al 

auggeated by the pulaar data, laply that either the denaltlea of the Ionized 

reglona are higher or the filling factor of Ionized gaa la larger than 

locally* Coablnlng thla Inforaatlon on ^n^> from pulaara with that 
2 

on <n > from the other observatlona which probe the Interatellar electron 
denalty should give a more complete and detailed picture of the structure and 
sourcea of Ionized gas .n the Galaxy. 
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TABLE 1 



120 pc. 


TABLE 2 


lA 

M 



(X^) “ probability that x* than the calculated 

value for n "degrees of freedou" 
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Figure It 


Figure 2t 


Figure 3: 


Figure 4: 


Figure 5: 


riCUII CAPTIONS 

Vertical dlaperelon Mesure (OH elnb) plotted egeinat horlaontal 
dleparjlon aaaaura (DM eosb) for a aaaple of 149 pulaara In the 
Inner Galaxy (310** < t < 50*^). Tha data la froa Manehaater et. 
el. 197P, Newton et> al. 1981 and Aahworth and Lyna 1981. 

Dlatrlbutiona of DM elnb for Molonglo aurvay pulaara with DM coab 
< 100 pc. ca~^ and DM coab > 100 pc. ca**^ for longitude rangee 
310° < t < 50° and 330° < 1 < 30°. 

Monte Carlo generated dlatrlbutiona of DM elnb for tha beet 
fitting electron model [cf. Eqn. (7) and Table 1]. 

Radial dependence of the dlak electron denalty, n^ (R), and of 
the central denalty, n^(R), and ecale height, h^(R), of the 
extended electron component for the beat fitting model. 

Monte Carlo generated dlatrlbutiona of DM elnb for an electron 
model with no radial dependence of the extended component [cf. 
Eqn. (9)]. 
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